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' Abstract. Using a combination of theory and computer simulations, we study the 

C/2 . translocation of an RNA molecule, pulled through a solid-state nanopore by an optical 

tweezer, as a method to determine its secondary structure. The resolution with which 
the elements of the secondary structure can be determined is limited by thermal 
fluctuations. We present a detailed study of these thermal fluctuations, including the 
' O ' frequency spectrum, and show that these rule out singlc-nucleotide resolution under 

^ . the experimental conditions which we simulated. Two possible ways to improve this 

Q I resolution are strong stretching of the RNA with a back-pulling voltage across the 

membrane, and stiffening of the translocated part of the RNA by biochemical means. 
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1. Introduction 

New developments in design and fabrication of nanometer-sized pores and etching 
methods, in recent times, have put translocation at the forefront of single-molecule 
experiments [1-9], with the hope that translocation may lead to cheaper and faster 
technology for the analysis of biomolecules. The underlying principle is that of a Coulter 
counter [10]: molecules suspended in an electrolyte solution pass through a narrow pore 
in a membrane. The electrical impedance of the pore increases with the entrance of a 
molecule as it displaces its own volume of solution. By applying a voltage across the 
pore, the passing molecules are detected as current dips. For nanometer-sized pores 
(slightly larger than the molecule's cross-section) the magnitude and the duration of 
these dips have proved to be rather effective in determining the size and length of these 
molecules [11]. 



Figure 1. The experiment in scliematics, illustrated by only 6 CG-bonds for clarity. 
An RNA molecule composed of N nucleotides is pulled through a solid-state nanopore 
in a membrane (placed at z — 0) towards the right {z > 0) using an optical tweezer, 
represented by a parabolic potential, at a constant speed vtw The bottom of the 
potential is located at ztw', and the latex bead is located at Zb- The number of 
monomers located on the trans-side of the membrane is called s. The monomers are 
numbered, starting from the end which is attached to the latex bead; consequently, the 
nucleotide located in the pore is labelled s. A potential difference 2V is also applied 
across the membrane. 



With a membrane placed at 2; = 0, a translocation experiment for determining the 
secondary structure of an RNA molecule [12-17,20] composed of N nucleotides proceeds 
as follows (see Fig. [1]). One end of the folded RNA, which is almost completely located 
on the left (cis) side {z < 0) of the membrane, is pulled through a solid-state nanopore 
to the right (trans) side [z > 0), and a latex bead is attached to nucleotide 1. An 
optical tweezer captures the bead, and pulls the RNA through the pore at a constant 
speed ftw The monomer number at any given time located inside the pore is denoted 
by s. We denote the voltages on the cis and trans side by +V and —V respectively; 
a potential difference 2V is thus applied across the membrane, to increase the tension 
in the translocated part of the RNA so that no secondary structure can form between 
the tweezer and the pore. During this process, the force exerted by the tweezer on the 
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RNA is monitored. Since the pore is narrow, for translocation to proceed, the bonds 
between the basepairs forming the secondary structures must be broken at the pore. The 
breaking of the basepair-bonds is detected as increased force on the tweezer. The force 
on the tweezer as a function of time can then be translated into the binding energies of 
the basepairs as a function of distance along the RNA, yielding a wealth of information 
on the secondary structure of the RNA. 

Note that actual RNA chains may be pulled through the pore either from the 3' or 
from the 5' end, and the respective force extension curves may pick up both the initial 
and final location of the stems. Prom both the force-extension curves, with the help of 
a probabilistic sequence alignment algorithm [15], one can subsequently reconstruct the 
base-pairing pattern, the success of which clearly depends on the accuracy with which 
the arrests in the force-extension curves — as experienced by the tweezer bead — can 
be tied to the actual translocation coordinates along the backbone of the RNA. 

Given that the distances are of nanometer-scale in this experiment, thermal 
fluctuations of the polymer are expected to blur the arrests of the force extension curves; 
i.e., blur the coherence between the force exerted by the tweezer and the nucleotide 
number located in the pore. Since in this experiment one cannot track the events at the 
pore, the unpredictability of the amount of low-frequency noise in solid state nanopores 
seem to be the main barrier to progress in this field [18, 19]. The central question 
addressed in this paper, therefore, is the level of resolution (in units of a nucleotide) 
that can be achieved by this experiment. All throughout this paper, wc define resolution 
as the accuracy with which the location of secondary structure can be determined along 
the backbone of the RNA — which is strongly affected by the coherence between the 
force exerted by the tweezer and the nucleotide number located in the pore. We address 
this question by studying the amplitude and frequency spectra of the fluctuations at 
the pore with a combination of theory and computer simulations. The amplitude 
spectrum determines the resolution limit that can be achieved by ensemble averaging, 
while frequency spectrum determines the resolution limit that can be achieved by time 
averaging. Note that the highest resolution that this setup can achieve depends on 
which of these two limits is higher; this prompts the study of both the amplitude and the 
frequency spectra of the fluctuations at the pore. Our study rules out single-nucleotide 
resolution under the experimental conditions which we simulated. Two possible ways 
to improve this resolution are strong stretching of the RNA with a back-pulling voltage 
across the membrane, and stiffening of the translocated part of the RNA by biochemical 
means. 

A related problem was considered by Thompson and Siggia [20], who studied 
whether a measurable signal can be obtained by pulling apart a DNA or RNA molecule 
by an atomic-force microscope. They formulated their theoretical analysis using an 
(equilibrium) partition sum that involved the interaction energy between the unzipped 
strands. In the case of pulling apart a molecule by translocation, the unzipped strands 
of the molecule are separated by an impenetrable membrane, so the unzipped strands 
cannot interact directly unless one of the strands translocates through the pore; i.e.. 
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the process of unzipping cannot be decoupled from the dynamics of translocation. 
Consequently, the method of Thompson and Siggia cannot be easily imported to 
study our setup. We also note that recently a number of researchers, e.g., Sauer- 
Budge et al. [21] and Bockelmann et al. [22] have studied the case of pulling apart 
a molecule by translocation: however, their formulations do not take into account any 
dynamics of translocation, and therefore they only provide a simplified analysis of the 
problem. Given that the (anomalous) dynamics of translocation involves long memory 
effects [27-31], we follow a different method in this paper; this allows us to study the full 
dynamic problem (i.e., including the frequency spectrum of the thermal fluctuations). 

Our work is related to the study of Bundschuh et al [16] of translocation of RNA 
or DNA through a nanopore, describing slow and fast regimes of translocation: for the 
former, the cis side of the RNA molecule essentially remains equilibrated at almost all 
times, while for the latter, the base-pairing pattern on the cis side is essentially frozen 
during unzipping. Our analysis describes a maximum pulling velocity of the optical 
tweezer that allows the trans side of the molecule enough time to always remain in its 
steady state, thereby providing the quantitative distinguishing characteristics between 
the two regimes described by Bundschuh et al . 

The structure of this paper is as follows. In Sec. [2], we describe our computer 
model. In Sees. [3] and HI we analyse the problem without and with thermal fluctuations 
respectively. In Sec. Owe conclude the paper with a discussion on the results. 

2. Computer model 

We model the RNA with N nucleotides as a lattice polymer with monomers on a face- 
centred-cubic lattice. Multiple occupation of the same site is forbidden, i.e., the polymer 
is described by a self-avoiding-walk. For practical purposes, this restriction is lifted for 
consecutive nucleotides along the chain. The dynamics of the polymer consists of single- 
nucleotide hops to nearest-neighbour sites, attempted at random with rate unity, and 
accepted with Metropolis probabilities. This model [23-26] describes both reptation and 
Rouse moves, but does not include explicit hydrodynamics. We have used this model 
successfully to simulate polymer translocation under various circumstances [27-31]. 
Since our model is a variant of a freely-jointed-chain, we expect it to reproduce poly(U) 
RNA behaviour reasonably well [32]. 

In translocation experiments with biological nanopores, e.g., alpha-haemolysin, the 
polymer might show sequence-specific binding and unbinding to the pore wall [6,7,33]. 
Such interactions between the polymer and the membrane are not expected to play 
a role in experiments of translocation through a solid-state nanopore, as used in the 
experiments of Refs. [8,9]. It has also been suggested that the translocation of single- 
stranded DNA through alpha-haemolysin nanopore is direction specific [34] (3' to 5' 
as opposed to 5' to 3'); in the same paper, by computer simulations, the authors 
demonstrated that such direction-specificity should not be present when the pore 
diameter is > 1.5 nm. Given that the typical diameters for solid-state nanopores are 
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> 2 nm [35], in our simulations, we neglect interactions between the polymer and the 
membrane, other than excluded- volume interactions (i.e., the polymer cannot cross the 
membrane other than through the pore). 

Since we want to study how secondary structure influences the translocation 
process, we add the ability for parts of the polymer to form two hairpins. The real 
RNA sequence which comes closest to our approach is a poly(U) RNA with a sequence 
composition f/30 {UqoGs2UqCs2)2 Uqo, wherein each nucleotide corresponds to a monomer: 
two C and G-nucleotides on neighbouring lattice sites can form a bond with an affinity 
EcG = 2.3 k^T, but we do not allow GU pairing. The latter is a simplification from how 
a real RNA molecule with the above sequence would behave, but with this simplification 
we a priori know what to expect for the secondary structure of this polymer — namely 
two hairpins, each with 32 CG-bonds [36] — on which we study the effect of thermal 
noise that limit the achievable basepair resolution of the secondary structures. 

Our model of the optical tweezer is that the latex bead, i.e., nucleotide 1 feels 
a spherically symmetric harmonic trap with spring constant ktw, centred around the 
location of the optical tweezer at a distance ^tw from the membrane. 

It is clear that our model does not capture the full details of a real laboratory 
experiment: indeed a more detailed model could include explicit hydrodynamics, 
detailed RNA interactions such as GU pairing, and a more elaborate description of 
the charge distributions on the RNA. Leaving out explicit hydrodynamics does alter 
polymeric motion (from Zimm to Rouse dynamics) and therefore is likely to affect the 
polymer's memory effects. Although at this moment we do not precisely know how the 
memory effects that are relevant for the present problem — discussed in Sec. H] — will 
alter when explicit hydrodynamics are incorporated, the low-frequency domination of 
the memory effects will not disappear. 

To correspond to experimental parameters we use a lattice spacing of A = 0.5 nm, 
comparable to the persistence length (as well as the typical inter-nucleotide distance) 
for poly(U) [37]. The resulting forces measured at the tweezer are 60 pN or less, similar 
to experimental values [37]. Equating the diffusion coefficient 2-5 x 10~^ cm^/s for 
Uq [38] to that of a polymer of length = 6 in our model, one unit of time in our 
simulations corresponds roughly to 100 ps. The tweezer velocity is one lattice spacing A 
per 300, 000 units of time, or ~ 20/xm/s, comparable to typical experimental velocities. 
The time scale A/wtw is larger than that (or of the order) of the longest time-scale of the 
translocated part of the RNA, implying that the translocated part of the RNA can be 
treated as properly thermalised at all times. In our simulations, the value of AU = 2qV 
ranges from 0.4 to 2.75 k-QT. Given that at room temperature k-^T = 25 meV and 
assuming that each nucleotide carries an effective charge around q = 0.5 times the 
electron charge (due to Manning condensation, which limits the charge to one electron 
charge per Bjerrum lengtljjf) [40], our simulations correspond to an experimentally 

f Due to Manning condensation, the effective charge per unit length is hmited to approximately one 
electron charge per Bjerrum length. In water, the Bjerrum length is about 7 A. Since the typical RNA 
base pair distance is « 3.4A, the effective charge is approximately 0.5e per nucleotide. In the pore, the 
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applied voltage differences ranging from 10 mV < < 70 mV. 

A typical simulation output is presented in Fig. [2l It consists of the force Ft^{t) 
exerted by the optical tweezer as a function of time: this information is readily accessible 
in real experiments. In our simulations, we also monitor the number s of monomer 
located in the pore as a function of time: this information is typically not accessible in 
real experiments. In the simulation of figure [2], the force exerted by the tweezer hovers 
around a fixed strength 3.4/cbT/A (approximately 30 pN), except between t = 1.4 x 10^ 
to 3.4 X 10^ resp. t = 4.6 x 10^ to 6.7 x 10^, when the first and second hairpins are 
pulled through the pore. Indeed, the top panel of Fig. [2] shows that at the onset of 
these intervals, s{t) is almost constant, around 90 and 220, the starting locations of the 
hairpins. 



2e7 4e7 6e7 




Time 



Figure 2. Upper panel: nucleotide in the pore s as a function of time for a poly(U) 
RNA of composition U^q (C/6oG'32C/6C'32)2 f^eoj pulled with constant velocity Vtw = A 
per 300,000 time steps (approximately 20 /im/s). The binding energy of each bond is 
set to EcG = 2.3 k^T, with 2qV = 1.5 /cbT'. Every data point is an average of 1500 
consecutive measurements each 100 time steps apart, with the standard deviation 
represented by the error bars. The two straight lines are guides for the eye. Lower 
panel: the corresponding chain tension measured by means of the optical tweezer, with 

dielectric constant is significantly lower than that of water, consequently the Bjerrum length (which 
is inversely proportional to the dielectric constant) is much larger, and hence the effective charge is 
much lower (a tenth of an electron charge per nucleotide or even less [21,39]). For our work, a key 
quantity is the stretching force, determined by the energy difference across the pore, set by the effective 
charge in solution. Thus, for our purpose, the relevant effective charge is 0.5e per nucleotide. The main 
consequence of the much lower effective charge inside the pore is that RNA is less eager to enter the 
pore, but that has no consequences for this work. 
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3. Translocation without thermal fluctuations 

First we discuss what sort of information on the secondary structure can ideally — i.e., 
in the absence of thermal fluctuations — be obtained from Ft„(t). We do this under the 
assumptions that the force extension curve of the RNA without any secondary structure 
is sequence-independent, and that the tweezer velocity is low enough for the force exerted 
by the tweezer to maintain a uniform chain tension all along the translocated part of 
the RNA. Then, (p is uniquely determined by the relative extension x = z^/s, and is 
balanced by Ft^it), i.e., 

4> = J^{x) = Ft^ = kt^{zt^ - Z],), (1) 

where fctw is the stiffness of the optical tweezer. Additionally, the equality of the rate of 
work done by the tweezer and the gain in free energy by the translocating nucleotides 
at the pore yields 

F^„dzy, = {AU -TAS)ds, (2) 

where AS* is the entropic cost per nucleotide translocation due to the imbalance of the 
(entropic) chain tension across the pore, and AU is the energetic cost per nucleotide 
translocation. If translocation of the nucleotides does not involve breaking of CG-bonds 
at the pore, then AU = AUc = ^qV, otherwise AU = AUb = 2qV + Ecg- Thus, given 
and {AU — TAS), Eqs. ([1]) and ([2]) determine both the tweezer force and the relative 
extension. 

During the translocation of the first 90 nucleotides of Uso {UQoG32UeC32)2 Uqo, no 
secondary structure is broken at the pore — consequently, AU = AUc — and the tweezer 
force remains constant at Ft^{t) = 3.4/cbT/A (approximately 30 pN). The speed of 
translocation is then given solely by Eq. ([1]), with ib = "^tw, as s = ftw/-^~M-^tw(^)] = 1 
nucleotide per 252, 000 time steps, for which we have used the numerically obtained 
force extension curve (inset. Fig. [3]); this speed s is shown in Fig. [2] by the upper (blue) 
line. 

Following the arrival of the first hairpin at the pore, translocation requires breaking 
of the CG-bonds, and consequently, AU increases by Eqg = 2.3 k-^T. Both the tweezer 
force and the relative extension adjust to new values, determined by Eqs. ([T]) and 
([2]) for the new value for AU. The resulting tweezer force equals Ft^(t) = l.bk-oT/X 
(approximately 70 pN); the correspondingly adjusted speed s is shown in Fig. [2] by the 
lower (red) line. 

After the translocation of the first 32 G-nucleotides, AU returns to its base value 
AUc- The tweezer force and the relative extension, too, fall back to their pre-hairpin 
values. This is seen in Fig. |2]by s{t) leaving the lower red line sharply to re-coincide with 
the upper blue line; i.e., quite a few nucleotides at the end of the hairpin translocate 
nearly immediately. This chain of events repeats itself during translocation of the 
second hairpin: first the tweezer force increases gradually to its higher value and the 
translocated distance approaches the lower red line, then the tweezer force decreases 
steeply to its lower value and the translocated distance jumps to the upper blue line. 
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In conclusion, most features of Fig. [2] are qualitatively well-understood. The above 
framework can be easily extended to a wider set of bond strengths and more elaborate 
secondary structures. Without thermal fluctuations, the setup is perfectly suited to 
determine the secondary structure up to the nucleotide resolution, under the restriction 
that the consecutive bonds along the backbone of the RNA are of increasing strength; 
if strong bonds are followed by weaker bonds that are not strong enough to halt the 
translocation process, the breaking of the weaker bonds will not be accompanied by 
an increase of and the experiment will reveal little information about these weaker 
bonds. 



4. Translocation with thermal fluctuations 



In reality, thermal fluctuations are omnipresent in this nano-scale experiment, and as 
argued earlier, for solid-state nanopores they are the dominant source of noise at the 
pore. In fact, it is precisely the (thermal) fluctuations in s(t) that serve to blur the 
coherence between -Ftw(t) and s{t), and thereby limit the resolution that can be achieved 
by this experiment. We will now study both the amplitude and the frequency spectrum 
of these thermal fluctuations in s{t). 
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Figure 3. u^, the mean square displacement of s vs. reduced extension x — z\j/ at 
constant ztw — 300 A. The chain tension (f) is slowly increased by changing AJ7 — 2qV 
from 0.40 to 2.75 /cbT. Each data point required 80 independent polymers and 
simulation times of 20 million time steps, with a measurement every 100 time steps. 
Data points from direct simulations and Eq. ^ are represented as black diamonds and 
red squares. The dashed lines are cubic splines. The error bars represent statistical 
errors only. Inset: Rescaled force-extension curve for our model. 



The amplitude asit, AU) of the fluctuations in s{t) is that of an entropic spring at 
fixed extension zi, with one end tethered at the tweezer, while the number of nucleotides 
in the spring are allowed to fluctuate through the pore. Now consider a different 
problem — an entropic spring with an average, but fluctuating extension Zb- From 
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the equipartition theorem, these fluctuations are given by {SzD = k-oT/cz^ [with spring 
constant = {dJ-'/dz]^)s = J-''{x)/s, in which x = z^/s is the relative extension]. For 
the present problem, such fluctuations in zi, can be thought of to be mediated by the 
fluctuations in s, yielding [41] 

2 

= s{t)[x^J^'{x)]'^ ksT. (3) 

X 

In Fig. [31 Eq. ([3]) is compared to the simulation results for several values of 
AU, with constant Ztw = 300 A. Note that AU only serves to set the value of 0. Of 
practical importance is the observation that according to Eq. ([3D the amplitude of 
the thermal fluctuations decreases with increasing tension. Thus, an increase of the 
applied voltage difference will reduce the thermal fluctuations in s thereby increasing 
the resolution of the secondary structure determination. The same effect can also be 
achieved by increasing the chain stiffness, thereby increasing jF'(x). In practice, this 
could for instance be realised through chemical means. For instance, it is known that 
the salinity affects the chain stiffness [37]. Also, certain proteins (such as RecA for 
ssDNA) can be added on the trans-side of the membrane, to increase the chain stiffness 
dramatically, while at the same time reducing secondary structure formation on the side 
of the tweezer. 



a',it,AU) = {Szl) 



ds 
dzu 




Figure 4. (Rescaled) power spectrum of s {t), Ss (/), versus (rescaled) frequency for 
AU = 1.5 ksT, and ztw/A=60, 90, 120, 150, 300, 450, 600, 750. Each curve is composed 
of statistics from 80 polymers for 40 million time steps, the value of s corresponding 
to each curve are shown in the Figure. The solid line /"'^/^ is added as a guide to 
the eye. The unit of s^f along the horizontal axis, 1/100 time steps, is approximately 
equal to 100 MHz. 

For the frequency spectrum of s(t), we return to Ref. [30]. Therein we showed that 
s{t) and the chain tension imbalance across the pore are related via a time-dependent 
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memory kernel a{t). This result, adapted to the notations in this paper, is given by 

s{t)= f rft'a(t-O[0(O-0.<o(t')], (4) 



where (l)z<o is the chain tension of the RNA at the z < side of the pore, and 
a(t) ~ exp[— t/rjr]; thus, the memory kernel a(t) is described by a power-law decay 
with exponent —3/2, up to some cut-off time rjr; this cut-off time is found to increase 
quadratically with the length of the translocated RNA: tjt ~ [30]. The immediate 
consequence of Eq. (jlj) is that 

m-Kt')) ~ it-tr'/'exp[-|t-t'i/^^], (5) 

following the fluctuation-dissipation theorem [42], i.e., 

{S{t)s{t')) = [ dti [ dt2{s{t^)s{t2)) 




\t — t \ ^ Tj: 

with some constant C . This autocorrelation function, for |t — t'| ^ r^r, has the same 
form as that of a fractional Brownian motion (FBM) [43,44] 

(s(t)s(t')) = Ci (|t|2^ + |tf ^ -\t- ^) , (7) 

with Hurst parameter Ji = 1/4, and some constant C\. Since the r.h.s. of Eq. ([7]) is not a 
function of it — t') alone, the spectral density of the FBM is not well-defined. However, 
by applying generalised concepts such as the Wigner-Ville spectrum [43], a limiting 
power spectrum can be obtained as follows: Ss{f) = with a = 2H + 1 = 3/2, with 
some constant C2. We therefore expect a power spectrum, with ~ l/rjc- ~ 1/s^, 

^'^^^ (I/I < /c) • 

As shown in Fig. HI after having rescaled Ss{f) and / by and respectively, 
the collapse of the power spectra confirms Eq. ([H]). The spectrum is dominated by the 
low frequency fluctuations of s(t), which are not easily suppressed by time averaging. 

The translocated part of the RNA can be considered to be in a steady state as 
long as s{t) < T^^, i.e., the time between the passage of consecutive nucleotides through 
the pore does not exceed the time for the longest memory in the translocated part of 
the RNA. This time for the longest memory is set by the longest length Smax of the 
translocated part of the RNA. If the tweezer is pulled with a velocity exceeding 

vt„.it, At/) = X ■ /, ~ ^'"^.f !fL (9) 

■^maxl'^B-' ) 

where Parseval's theorem has been used to equate Eq. ([8]) to Eq. ([3]), the trans side 
of the molecule is not allowed enough time to reach the steady state. This results 
in forces at the tweezer that far exceed the forces required to pull monomers through 
the pore at steady state (of the molecule on the trans side). Fig. 4 shows that if a 
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potential difference of AU = l.^k-^T is applied, the power rescaled spectrum crosses 
over from a frequency-independent regime to a power-law regime if s^/ = 0.1 to 0.15; 
with a pulling velocity Vtw = A/300,000 (roughly ~ 20/im/s) the typical frequency is 
/ ~ 300, 000~^ ~ 0.3 MHz; the rescaled quantity s^/ then reaches the crossover once the 
translocated part of the chain has a length Smax ~ 200 nucleotides. Beyond this length, 
the translocated part of the chain will be pulled out of its equilibrium shape, and the 
force felt by the tweezer is no longer determined by what happens near the pore, but 
rather by what happens in the translocated part of the chain, with the consequence that 
limited information about possible secondary structure will be obtained. In targeted 
simulations we have indeed observed the behaviour that without any base-pairing, for 
tweezer velocities higher than A/300,000 (^20 fim/s), the force at the tweezer is an 
increasing function of time significantly once Smax ~ 200 has been reached. 

5. Discussion and conclusions 

Having demonstrated the importance of thermal fluctuations in s{t), let us now revisit 
the translocation of U30 {UqqG 32UgC 32)2^^0, cf. Fig. [2j Except between times 1.4 
and 3.4 x 10^ (resp. 4.6 and 6.7 x 10^), Ftw(t) is roughly SAksT/X (approximately 
30 pN), with fluctuations {AF^^(t)) ~ k^^k^T, while Eq. ([3]) explains the growth of 
(j1{t) linearly with s{t) (hence linearly with t). It is important to stress that while 
the fluctuations in Ftw(t) may be removed to a reasonable degree by time averaging 
(e.g., by reducing the velocity of the tweezer), the fluctuations in s(t) fundamentally 
reduce the accuracy with which RNA secondary structure can be observed. As 
soon as the hairpin is within a distance of order as{t,AU) from the pore, i.e., 
(sit)) + as(t,AUc) = 90 (resp. 220), the tweezer force increases already, to reach a 
plateau once {s(t)) — as(t,AUb) = 90 (resp. 220). Translocation at the pore then 
proceeds with a constant rate ftw/^~"'^(-^tw)5 until {sit)) + asit,AUb) = 122 (resp. 252) 
(in which we assume that the effective affinity Ecg is not altered much as the 
unzipping of the hairpin proceeds), at which time the chain tension is released by 
rapid translocation of the remaining nucleotides. These considerations clearly establish 
that only secondary structure elements with at least two times as{t,AUb) basepairs 
can be detected accurately, assuming that the affinities of consecutive basepairs are 
more or less alike. With the choice of parameters in the simulations leading to Fig. [2] 
which are comparable to typical experimental values, this means that the resolution is 
limited to ~ 8 nucleotides. Information on smaller lengths is washed out by the thermal 
fluctuations, making it very hard to retrieve. Moreover, the basepairs of the preceding 
stem of the secondary structure must be characterised either by weaker affinities, or by 
strong heterogeneity. 

Although the fluctuations in s(t) are not easily suppressed, Eq. ([3]) does leave open 
avenues for higher accuracy, either by increasing the relative extension x = Zb/ s, or by 
increasing the stiffness J^'{x) of the polymer on the trans side of the membrane. The 
former can be achieved by applying a stronger potential difference 2V (the strength 
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of which is of course hmited by experimental considerations). On the other hand, 
the polymer's stiffness can be actively enhanced, e.g., by altering the salinity of the 
solution [37]; alternatively, if one is interested in doing single-molecule experiments 
with ssDNA, the addition of RecA proteins to the solution only on the trans of the 
membrane may be of help as well (this will increase the stiffness of the polymer on the 
trans side of the membrane, but since these proteins cannot pass through the pore, the 
secondary structure on the cis side will be left unaffected) . 

If biological pores such as alpha-haemolysin are used for this experiment instead 
of solid-state nanopores, the translocating polymer will show sequence-dependent 
binding/unbinding to the pore wall. Although we cannot oversee all consequences of 
such interactions, we do not expect these to affect the coherence between the force felt at 
the tweezer and the translocated length, nor the frequency spectrum of the fluctuations 
in the translocated part of the chain. Hence, our conclusions should also apply to 
biological pores. 
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